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ABSTRACT 

We used the Project 1640 near-infrared coronagraph and integral field spectro- 
graph to observe 19 young solar type stars. Five of these stars are known binary 
stars and we detected the late-type secondaries and were able to measure their 
JH spectra with a resolution of R ~ 30. The reduced, extracted, and calibrated 
spectra were compared to template spectra from the IRTF spectral library. With 
this comparison we test the accuracy and consistency of spectral type determina- 
tion with the low-resolution near-infrared spectra from P1640. Additionally, we 
determine effective temperature and surface gravity of the companions by fitting 
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synthetic spectra calculated with the PHOENIX model atmosphere code. We also 
present several new epochs of astrometry of each of the systems. Together these 
data increase our knowledge and understanding of the stellar make up of these 
systems. In addition to the astronomical results, the analysis presented helps 
validate the Project 1640 data reduction and spectral extraction processes and 
the utility of low-resolution, near-infrared spectra for characterizing late-type 
companions in multiple systems. 

Subject headings: binaries: visual - instrumentation: adaptive optics - stars: 
individual (HD 77407, HD 91782, HD 112196, HD 129333, HD 135363) stars: 
late-type 



Introduction 



Project 1640 is a near-infrared integral field spectrograph (IFS) behind an apodized pupil 
Lyot coronagraph at Palomar Observatory. The instrument is fully detailed in lHinkley et al. 
( 1201 11 ). Between 2008 March and 2010 July, we carried out the Phase I survey, coupling 
the Project 1640 instrument with th e 241-element PALAO adaptive optics (AO) system 
( iDekany et al.Hl 9971 ; iTroy et al.l 120001 ) on the 5.1 m Hale Telescope at Palomar Observatory. 
For the Phase II sury ey, the instrument is mated with the high-order AO system PALM-3000 
( IBouchez et al.ll2010l ) on the same telescope; the Phase II survey started in 2011 October. In 
combination with a post-coronagraph ic wavefront calibration unit and image post processing 
( ICrepp et al.ll2011uPueyo et aU201l[ ). we expect to generate contrasts of 10~ 7 at a separation 
of 1". The project's goal is the detection and characterization of exoplanets and brown 
dwarfs. 

The stars in this paper were targets of the Formation and Ev olution of Planetary Systems 
(FEPS) Legacy Science program for the Spitzer space telescope (IMeyer et al.ll2004l ). selected 
to be solar-type stars spanning a range of ages from ~3 Myr t o ~3 Gyr. Those accessible from 
the northern hemisphere were surveyed for multiplicity by iMetchev fc Hillenbrand! ( 120091 ) . 
During the P1640 Phase I survey, we observed a number of young FEPS targets, as their 
youth makes them promising candidates for searches for exoplanets. This paper discusses 
five stars ((HD 77407, HD 91782, HD 112196, HD 129333, HD 135363) that have known 
stellar companions detected by previous AO observations. It also briefly discusses the stars 
that we did not detect any companions to. 



This paper furthers our understanding of these five stellar systems by adding to the 
known astrometry, which will eventually allow for the computation of an orbital solution. 
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We also provide estimates of the spectral types of the companions. Coupled with the even- 
tual mass estimates of the various spectral types, this will provide a measurement of the 
mass ratio and contribute to our understanding of the mass ratio distribution for solar type 
stars. The companion mass distribution is a key test of various stellar formation models 
( IKouwenhoven et all 12003 ; iReggiani fc Meverll201lh . 



In Section [2] we describe how the observations were carried out. We describe how the 
data were reduced and how the astrometry and spectrum of the companion were produced 
in Section El Section HJ discusses the results for each systems. The results of a model 
based method for determining the effective temperatures and gravity of the companions 
are presented in Section |5j The observations of stars that did not detect a companion are 
discussed in Section [6J Finally, Section [7] gives our summary and conclusions. 



New Observations and Previous Data 



All the stars were observed during the same observing run in 2009 March. HD 129333 
was observed on 2009 March 14 UT; HD 91782 and HD 135363 were observed on 2009 March 
16 UT. Finally, HD 77407 and HD 112196 were observed on 2009 March 17 UT. The mean 
seeing for each night was visually estimated to be approximately 1", which is slightly better 
than the median seeing for the observatory. The data collection plan was the same for all 
stars. First, several unocculted IFS images were acquired. These were mainly used to guide 
the star onto the occulting spot. After the star is occulted, multiple images were collected, 
each with an exposure time of 127.5 s. We collected 15 images of HD 77407, HD 91782 
and HD 112196. For HD 129333 we collected 10 image and for HD 135363 we collected 
only three images . Afte r collection, the data were processed with the method detailed in 
Zimmerman et al.l ( 1201 ll ). The processing turns each IFS image into a three-dimensional 
data cube, with 23 narrow-band images with central wavelengths between 1.100/xm and 
1.760/mi at 30 nm increments. 

We combine our P1640 astrometric measurements with multi-epoch Palomar/PHARO 
and K eck/NIRC2 astrometric data collected between 2002-2004 from the survey of iMetchev 
( 120061 ) . These data were taken with the Palomar 5m telescope usi ng PALAO with the 
Palomar High- Angular Resolution Qbser yer (PHARO) NIR cam era (IHayward et al.l 120011 ) 
and with the Keck II facility AO system ( jWizinowich et al.ll2000l ) using the NIRC2 camera. 
The PHARO distortion and pixel scale were calibrated to a precision of 0.15% t hrough a 
combi nation of astrometric mask and binary star measurements, as described in IMetchev 
( 120061 ). Additional information on the observations and data reduction for these data is 



found in IMetchev fc Hillenbrand! (12009( 1 
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3. Data Analysis 



Extracting spectra from IFS instruments is non-trivial with many challenging aspects 
including image registrat ion and cross talk between channels. These have been described in 
Zimmerman et al.l (1201 if ) which discusses the current P1640 data pipeline. Pipeline process- 
ing produces an image of the object in each of 23 wavebands resulting in a data cube. For 
the unocculted data cubes, we measured the binary star astrometry on each of the individual 
image frames with the program FITSTARS; it uses an iterative blind- deconvolution that fits 



the location of delta functions and their relative intensity to the data (jten Brummelaar et al. 



19961 . 12000J). After throwing out any image frames that failed to converge to a physical so- 
lution, the position angle and separation was computed as the weighted mean of the results 
from all the frames. The weights were set equal to the inverse of the RMS residual of the fit, 
a standard output of the FITSTARS program. The error bars for position angle and sepa- 
ration are the weighted standard deviation of the results. The plate scal e and the position 
angle offset were assumed to be constant between observing runs. See IZimmerman et al. 
( I2010l ) for more details on the plate scale of the Project 1640 instrument. 



The measured astrometry is shown in Table [TJ The table lists the star's Washington 
Double Star Catalog number, the HD number, the date of the observation in Besselian 
years, the position angle, separation and the instr ument that was used to collect the data. 
We also include the astrometry from the survey of iMetchev fc Hillenbrand! (120091 ) . In that 
paper, only the discovery astrometry was published; we include all the astrometry from the 
survey. For the PHARO data, we have excluded the 0.12 deg la uncertainty in the detector 
position angle, as that uncertainty is systematic and masks the already apparent orbital 
motions, especially in HP 112196 AB. The astrome try published here is slightly different 
than that published in IMetchev &: Hillenbrand! ( 120091 ) . In IMetchev! (120060 one of the corners 
of the PHARO detector was used as the origin for solving for the principal magnification 
component of the astrometric distortion. For the sake of absolute astrometric calibration, 
and for consistency with all other parts of the astrometric distortion analysis, the current 
data analysis used the center of the array as the origin of the distortion. This results in a 
small improvement in the astrometry. 

To create a spectrum of each object, we performed aperture photometry on each of the 
23 images in the data cube. This was done w ith the aver, p ro routine which is part of the IDL 
astroliqj and is an adaptation of DAOphot ( Stetson! Il9871 ). An aperture was drawn around 
each object as well as an annulus. The radius of the aperture and annulus was adjusted for 
each star, so it avoided the other star and avoided the occultation spot. The background was 



1 http://idlastro. gsfc.nasa.gov 
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set equal to the average intensity for each pixel in the annulus and was subtracted from each 
pixel in the aperture. The remaining power was then totaled. The total power as function of 
wavelength is the corresponding spectrum. The aper.pro code also produces an measurement 
error bar based on photon statistics and the error in the background estimation. 

The extracted spectrum is a convolution of the object spectrum, the instrumental 
spectral response function and the absorption due to the Earth's atmosphere. To remove 
the instrumental response and the atmospheric absorption, we used the primary in un- 
occulted images as a reference source. We used a spectrum corresponding to the spec- 
tral type of the primary from t he Infrared Telesco pe Facility (IRTF) Spectral Library 



( ICushing. Rayner fc Vaccal 120051 ; iRayner et al.l 120091 ). We divide the extracted spectrum 



of the primary by the reference spectrum; the result is the spectral response function (SRF) 
of the instrument and the atmospheric absorption. The companion's measured spectrum is 
then divided by the SRF and yields the calibrated spectrum of the companion. Even though 
the secondaries are visible in the unocculted images, we only measured their spectra in the 
occulted images, as it had much higher signal to noise. 

The error bars on the final reduced spectra are the combination of the errors of the mea- 
sured science spectra and the error on the SRF. The SRF errors come from the combination 
of the errors of the measured calibration spectra and the errors on the template spectra. The 
error bars of the measured spectra are dominated by the cr oss talk error, where ligh t from 



adjacent microlenses in the IFS overlap on the focal plane ( I Zimmerman et al.l 120 111 ). The 
measurement error from the aperture photometry is also very small. We had multiple data 
cubes of the occulted star, which reveal the companion. In these cases, the read noise was 
reduced by averaging the spectra from multiple data cubes. For the calibration observations, 
where we only had one or two observations, this was not the case and the error bars on the 
SRF are much higher. This turns out to be the dominant error term in the error bars of the 
final spectra. 

There is one astrophysical error term that is independent of our instrument. Using 
an IRTF standard star in the computation of the SRF rather than the true spectra of the 
primary causes a small error. There are two components of this error. One is that there may 
be an error in the reported spectral type of the primary; it may be either later or earlier than 
what we use. On most stars, the difference between various determinations of the spectral 
class is only a few subclasses. In these cases there is no change in the result. A larger error 
in the spectral type of the primary would result in a corresponding error in the secondary's 
spectral type. 

The other error component is that the template star will not have the exact same age 
or metallicity as the science target, and will have a slightly different spectrum. To estimate 
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the size of this error, we found two cases, where there were multiple observations of the same 
spectra type in the IRTF templates. We computed the ratio of the spectra from the same 
spectral type and plotted the result in Figure [TJ The error was ±4%, which is smaller than 
the cross talk error. Also the exact shape of the template error is unknown and will vary 
depending on close the parameters of the science target match the template star. For these 
reasons, it was left out of the formal error. Both components of the error can be eliminated 
if the spectrum of the primary is known from observations with other instruments. 



Our earlier papers iHinkley et al.l ( 120101 ) and IZimmerman et al.l ( 120101 ) did not calibrate 
the SRF on the target star, but instead used a separate observation of a calibration star to 
calibrate the SRF. As such, they had large errors where telluric water bands occur, between 
1.37-1.43 ^tm, and this portion of the spectrum was not used. Our hope was that the self 
calibration process using the method described above would improve the removal of telluric 
lines and allow for more of the spectrum to be used. This has proven to be moderately 
successful for three of the stars. For these stars, the spectra have increased error bars in 
the water ban. As stated above, the error bars are dominated by the SRF and the SRF is 
dominated by the signal to noise which is low in the water band. 

For two of the stars, self calibration did not work as well as the others. For HD 135363, 
the separation between the primary and secondary was small enough that there was spectral 
contamination, which affects calibration. In this case, the SRF had to be taken from another 
star. In HD 112196 the reason for the poor correction the spectrum of HD 112196 is unknown. 
All of the recorded spectra show a residual spectral absorption feature in the 1.37-1.43 
/iin water band, although only three minutes separate the unocculted and the first occulted 
image. So it is probably not a case where the humidity changed during the observations. In 
the end, we did not use the portion of the HD 112196 spectra in the water band. 

After the spectrum of the companion was extracted, it was then compared against the 
spectra in the IRTF Spectral Library, which were used as templates. Since the primaries 
are all young stars, we compared the spectrum against the FGKM main sequence stars and 
LT brown dwarfs. The template spectra were binned and smoothed in order to pr oduce 



the equivalent spectra to having the star observed by P1640 (IZimmerman et al.ll201lf ). The 
spectral standards were then compared against the measured spectrum. The metric for 
comparison was the weighted sum of the squares of the residual (SSE). 



SSE = J2MSx-Rx) 2 , (1) 

A 

where, S\ is the measured spectrum at a given wavelength, R\ is the binned reference 
spectrum at the same wavelength and w\ is the weight at the wavelength. The best fit 
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Fig. 1. — The solid line (orange in online edition) is the ratio of spectra for two F8V stars, 
and the dashed line (green in online edition) is the ratio of two G8V stars. All four spectra 
were taken from the IRTF Spectral Library. The ratio of the spectra is the relative error 
that would be caused by using one star's spectrum to calibrate the SRF of the other. All 
spectra were binned and smoothed to appear as they would if taken by P1640. 
(A color version of this figure is available in the online journal.) 

reference spectrum was the one with the minimum value to the metric. The weights were 
set equal to the inverse of the computed error at each wavelength point. Usually several 
adjacent spectral types fit equally well within the error bars and produce a range of spectral 
types. 



Results 



4.1. HD 77407 



HD 77407 (WDS 09035+3750) is a binary system of a young GOV and a late type 
companion. It is a nearby s tar, with a distance of 30.1 pc ( Ivan Leeuwenl 120071 ). The com- 
panion was first detected by iMugrauer et al.l (120041 ) via near-IR AO. They also measured its 
radial velocity over two years and show the star exhibits a long term radial velocity trend. 
The radial velocities have a great deal of noise, presumably from chromospheric activity 
of the young star. Due to this noise, they are able to place only the broadest constraints 
on the mass of the companion of 0.3-0.5 M Q . They estimate the companion spectral type 
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to be M0V-M3V based on ij- band photometry and evolutionary models of late-type stars. 



Metchev &: Hillenbrandl ( 120091 ) confirmed that the objects share common proper motion. 



The extracted spectrum of HD 77407 B is shown in Figured) overplotted are the spectra 
for M2V, M4V, M6V and M8V. These show how the spectra vary over over the range of 
spectral types. Unfortunately the companion landed near the edge of the PI 640 field of 
view. A portion of the third Airy ring is cut off in the longer wavelength images in the data 
cube. This adds additional uncertainty to the extracted spectra. The best fit (via Equation 
[T]) f or the companion's sp ectral type is M3V-M6V. This is later than the M0V-M3V derived 
by iMugrauer et al.l (120041 ) , but that determination was done with onlyif-band photometry. 
Still our results on HD 77407 B are the most uncertain of any of the stars in this paper. 
Additional observations with either higher resolution or at a longer wavelength would provide 
more information. 



The astrometric measurements from P1640 and from iMetchevi (12006!) ar e show n in 
Table [p They a r e cons istent with the prior measurements from IMugrauer et al. ( 2004 ) and 



Lafreniere et al.l ( 120071 ). Using stellar masses corresponding to the spectral types of the 



primary and the secondary, and assuming a circular face-on orbit we estimated the orbital 
period of ~300 years. This is an upper bound to the orbit and depending on the eccentricity 
of the orbit, could be significantly shorter. 



4.2. HD 91782 



HD 91782 (WDS 10368+4 743) appears to be a nondescript GO star with few published 
details. Its distance is 61.4 pc (Ivan Leeuwenl 120071 ). The FEPS project estimated the age to 
be 97 Myr by an average of thr ee techniques: rotation rate, ac tivity, and lithium abundance. 
A companion was detected by iMetchev fc Hillenbrandl (120091 ) and shown to have common 
proper motion. 

In the P1640 data, the companion to HD 91782 is much fainter than the primary, and it 
is difficult to detect the companion in many frames, especially around the water absorption 
bands. Only 13 images produced meaningful results and they have a large scatter in the 
measured astrometry resulting in the large error bars in Table [H The astrometry from 
Metchevi (120061 ) are also given in Table [TJ 



There is no published luminosity class for HD 91782 (jKharchenkol 1200 It IWhite et al. 



20071 ). Based on the age estimate, we used a main sequence GO star from the IRTF Spectral 
Library as the calibration spectrum. The extracted spectrum for the secondary is shown in 
Figure The best fit for the spectral type of the companion is M9V. 
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Fig. 2. — The data points are the measured spectrum of HD 77407 B. The black solid line is 
the spectrum for an M2V star, the dotted line (green in online version) is for an M4.0V star, 
the dashed line (red in online version) is of an M6.0V star and the dot-dashed line (orange 
in online version) is for a M8V star. 

(A color version of this figure is available in the online journal.) 

The astrometry for HD 91782 is shown in Table [0 Only a small fraction of the orbit 
has been observed. In order to get an upper bound to the orbital period, we assume it has 
a circular orbit with a inclination of zero. Using Kepler's third law, we get a period of ^450 
years. 



4.3. HD 112196 



H P 112196 (WPS 12 547+2206) is a F8V star (IHarlan fc Tavlorlll970h at a distance of 
34 pc (Ivan Leeuwen! 120071 ). The FEPS project estimated the age to be 134 Myr (an average 
of calcium HK and lithium abundance methods), though the ro tation-speed age measurement 
produ ced an anomalous age of 2.5 Gyr. Using isochrone fitting. iHolmberg. Nordstrom fc Andersen 
(120091 ) estimated the age to be 4.9 Gyr, with a la limit of 2.4 to 7.1 Gyr. 



A companion was detected by iMetchev fc Hillenbrand! (120091 ) . Using observations at 
multiple epochs they concluded it was a physical companion b ased on common p roper mo- 
tion. The star was observed by speckle interferometry in 1999 (IMason et al.ll200ll ) and 2001 
(Mason et al. in preparation), but the companion was not detected due to the faintness of 
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Fig. 3. — The data points are the measured spectrum of HD 91782 B. The black solid line is 
the spectrum for an M2V star, the dotted line (green in online version) is for an M4.0V star, 
the dashed line (red in online version) is of an M6.0V star and the dot-dashed line (orange 
in online version) is for a M8V star. 

(A color version of this figure is available in the online journal.) 



the companion. 

There was a single unocculted data cube of HD 112196 where both components can be 
detected. The companion is easily identified in each of the data cube slice. We measured 
the astrometry of in each channel of the data cube with FITSTARS as described in Section 
HI The data channels at wavelengths of 1100 nm and 1130 nm failed to converge. These two 
channels suffered from poor signal to noise due to the atmospheric water b and abso r ption . 
The resulting astrometry, as well as the astrometry from the observations in iMetchevI ( 2006 ) 
are shown in Table [TJ 

Unlike the other stars in this paper, the extracted spectrum for HD 112196 B shows 
signs of telluric water absorption after the SRF calibration process. The reason for this is 
unknown. We examined the spectrum from each of the individual data cubes and they were 
consistent, which seems to eliminate the possibility that the water column density changed 
during the observations. We did not use those data points in our spectra fitting process. 
The measured spectrum for HD 112196 B is in Figure HI The best fit for the spectral type 
of the companion is M2V-M3V. 



The astrometry for HD 112196 is shown in Table [TJ The upper bound estimate for the 
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orbital period yields a period of ^320 years. 
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Fig. 4. — The data points are the measured spectrum of HD 112196 B. The black solid line 
is the spectrum for an M2V star, the dotted line (green in online version) is for an M4.0V 
star, the dashed line (red in online version) is of an M6.0V star and the dot-dashed line 
(orange in online version) is for a M8V star. 
(A color version of this figure is available in the online journal.) 



4.4. HD 129333 



HD 129333 (WDS 14390+6417, EK Dra) is a young solar analogue (IDorren & Guinan 



1994 ) t hat is highly variable and is possibly a member of the Pl eiades moving group (ISoderblom &: Clements 
19871 ). iMontes et ajj ( 1200 ll ) gives a spectral type of G1.5V. iDuquennoy et al.l (Il99lh were 
the first to deter mine that this was a binary system via radial velocity measurements, and 
shortly thereafter iDuauennov &: Mavorl (ll99ll ) published a spectroscopic orbit of the system. 
Metchev fc Hillenbrand! (120041 ) first reported a direct image of the companion and confirmed 
its phy sical associat i on fro m AO observations with the Palomar Hale telescope. Indepen- 
dently, iKonig et al.l ( 120051 ) reported multi-epoch speckle interferometry observations of the 
companion, which they combined with radial velocity measurements to compute an orbital 
solution for the binary. They estimate the period to be 45± 5 years with an eccentricity of 
0.82±0.03. They derive a 0.5±0.1 M m ass for the secondary. Near-IR (J and K band) spec- 
troscopy of HD 129333 B published in iMetchev &: Hillenbrand! (120041 ) sets the companion 
spectral type at M2±1V. 
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The extracted P1640 spectrum of the companion is shown in Figure [51 The best fit for 
the spectral type of the companion is in the range M2V-M3.5V. This is in agreement with 
the M2V±1 spectral determination of iMetchev fc Hillenbrand! ( 120041 ). who produced their 
result using much higher resolution spectrum (R=1000). 
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Fig. 5. — The data points are our measured spectrum of HD 129333 B. The black solid line 
is the spectrum for an M2V star, the dotted line (green in online version) is for an M4.0V 
star, the dashed line (red in online version) is of an M6.0V star and the dot-dashed line 
(orange in online version) is for a M8V star. 
(A color version of this figure is available in the online journal.) 

The measured astr ometry is shown in Table [TJ The orbital residuals between our data 



points and the orbit of iKonig et al.l (120051 ) are shown in Table [2 The 2009 observation has 
the largest difference in position angle, but the smallest difference in separation. It is hard to 
say what fraction of the difference is due to errors in the computed orbit and which is due to 
errors in the P1640 astrometry. We will continue to monitor the astrometric measurements 
from P1640 to determine if there are any systematic errors. 



4.5. HD 135363 



HD 135363 ( WPS 15079+7612) is a nearby young G 5V star (IMontes et al.ll200lh . Its 
distance is 29 pc (Ivan Leeuwenl 120071 ). iBubar et al.l (120071 ) give an age of 36 lg 4 Myr and a 
mass of 0.78±0.01 M based on its location on the HR diagram. They determine th at it is 
a post T-Tauri star. It is a member of the IC 2391 supercluster (IMontes et al.l 1200 ll ). 
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The system has multiple components , of wh ich the physical nature is still uncertain for 
two of the pairs. iMetchev fc Hillenbrand! ( 120091 ) reso lved the Aa,Ab pair and sh owed that 
the two components shared common proper motion. iLepine fc Bongiornol ( 120071 ) identified 
a common proper motion companion with a separa tion of 116^5, th e B component. The C 
compo nent with a separation of 7'.'5, was detected by lLafreniere et al.l (120071 ) . lLowrance et al. 
( 120051 ) observed the star with NICMOS in 1998 and failed to resolve the Aa,Ab pair. The 
NICMOS coronagraph has an inner working angle of OM, and the companion was most likely 
within this angle. They detected another possible companion with a separation of 17'.'11, 
the D component. The C and D components have not been confirmed and the question 
remains if they are gravitationally bound to the primary. The star was observed by speckle 
interferometry in 2001 (Mason et al. in preparation), but no companions were detected due 
to either the contrast exceeding the dynamic range of the instrument or to the system having 
a separation wider than their speckle camera's 3" field of view. 

The SRF extracted from the HD 135363 images produced unphysical results. This is 
probably because of contamination from the secondary. Instead the SRF from HD 91782 
taken the same night was used. This produced the spectrum of the secondary shown in 
Figure |6J The best fit for the spectral type of the companion is in the range of M2V-M4V. 
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Fig. 6. — The data points are our measured spectrum of HD 135363 Ab. The black solid 
line is the spectrum for an M2V star, the dotted line (green in online version) is for an 
M4.0V star, the dashed line (red in online version) is of an M6.0V star and the dot-dashed 
line (orange in online version) is for a M8V star. 
(A color version of this figure is available in the online journal.) 
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The measured astrometry is shown in Table [TJ From Hipparcos astrometry. iMakarov fc Kaplan 



( 120051 ) derive a binary star period of 39.8 yrs, presumably this is from the Aa,Ab pair. Our 
measurements show the companion only moved 12° in seven years. An estimate of a circular 
face-on orbit gives a period of 29 yr. If it was such an orbit, we would have expected to 
see 24% of the orbit or 87° in the seven years of observations. Of the stars in this paper, 
HD 135363 is the highest priority for additional astrometry. The period is probably on the 
order of 40 years and seven years of that has already been observed. Observations every few 
years over the next few decades will enable an orbital computation. 



5. Model Comparisons 



In addition to the spectral type analysis described above, we attempt to infer the physical 
properties of the compa nions from best fit syn thetic spectra calculated with the PHOENIX 
model atmosphere code (IHauschildt et al.lll999l ). The atmosphere struc tures (temperature 
and pressure as a function of optical depth) were calculated as described in lRice et all ( l2010ai ) 
for a larger range of effective temperatures, in this case 1400 to 4500 K in 50 K increments. 
The range of surface gravities is 3.0 to 6.0 dex in cgs units), calculated in 0.1 dex increments. 
The PHOENIX code produces spectra with 4 A resolution at near-infrared wavelengths, which 
is more than adequate for the low-resolution spectra (8X ~300 A) obtained with Project 1640. 
These data constitute the pre-calculated grid of synthetic spectra referenced below. 



5.1. Fitting Procedure 



The fitting procedure we used is based on the methods described by lRice et al.l ( )2010af ) 
for high-resolution spectra of very low mass stars and brown dwarfs. In addition to the larger 
range of effective temperatures considered for this study, we have improved the sophistication 
and efficiency of the procedure and incorporated corrections for known systematic errors 
intro duced by the Project 1640 spectral extraction procedure (described in lZimmerman et al. 
20111 ). A more detailed description of the fitting procedure and the results of its application 
to spectral templates from late-type objects in the IRTF spectral library will be presented 
in Rice et al. (2012, in prep.). 

The best-fit parameters and their uncertainties are obtained in th r ee ma in steps. First, 
a goodness of fit parameter similar to that defined in ICushing et al.l (120081 ) is determined 
for every pre-calculated synthetic spectrum in the model grid, nearly 2,000 spectra in total. 
The goodness of fit values are used to constrain the range of temperatures and gravities for 
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the least-squared minimization calculations in the second part of the fitting procedure. The 
range of the constrained parameters varies for each object depending on the topology of the 
goodness-of-fit surface. 

The second part of the code introduces linear interpolation in flux between pre-calculated 
models, resulting in higher-precision best-fit parameter results than would be possible with 
only the pre-calculated model grid. The initial parameters of the least-squares minimizations 
are randomly assigned from the constrained range of temperatures and gravitie s. The least 



squar es minimization is accomplished in 10 5 iterations of the MPFIT code (IMarkwardt 



20091 ). resulting in a distribution of best fit parameters and corresponding % 2 - values. The 
parameters with the lowest \ 2 are adopted as the starting point for the third step. 

Finally, a Markov Chain Monte Carlo (MCMC) is conducted for the entire range of 
parameters in the model grid, and again the code linearly interpolates in flux in order 
to analyze an effectively continuous grid of models. The MCMC chain for each observed 
spectrum is 10 6 links, with the first 10% considered "burn-in" and not included in the final 
analysis. The "jump size" (width of the normal distribution from which the new parameter 
was randomly selected) was optimized to produce an acceptance rate (number of accepted 
jumps divided by the total number of links) of 0.3-0.4 and varied from 100 to 750 K and 0.2 
to 1.5 dex in surface gravity. The best-fit parameters are the modes of the 9xl0 5 post-burn 
links, which are calculated from the two-dimensional density distribution. The uncertainties 
are the l-a width of the density distribution in each parameter. 



5.2. Fitting Results 

Table [3] shows our estimates for T e fff, log(g), and uncertainties inferred from the model 
fitting as well as the spectral types we determined from IRTF spectra. The results support 
potential success of efficiently and accurately inferring the physical properties of late-type 
companions directly from low-resolution spectra. However, the results also reveal systematic 
errors that need to be examined before the procedure can be reliably applied to completely 
unknown companions. 

One object with excellent data, HD 129333 B, resulted in a best-fit temperature that is 
well-determined and consistent with the object's spectral type and age. The best-fit surface 
gravity is at the high end of the allowed range, and is unphysical according to evolutionary 
models. The results for this object are discussed further in Section 15.31 



For the other objects, best-fit temperatures tend to be higher than those inferred from 
spectral type-effective temperature conversions of late- type dwarf objects (e.g., (ILuhman et al 
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19991 . 120031 ; iGray &: Corballyl 120091 ) ). although the model fitting temperature uncertainties 
are as high as ±320 K. The best-fit temperatures of three objects (HD 77407 B, 129333 B, 
and HD 135363 Ab) are consistent with the spectral type determined from the P1640 data. 
HD 135363 Ab produced best-fit surface gravities of 5.90 dex, at the upper boundary of the 
models, and a slightly higher temperature than implied by its spectral type, although within 
the large uncertainty. HD 112196 B has the largest uncertainty in temperature of any object, 
and the MCMC distribution extends to ~3500 K, consistent with the earliest spectral type 
estimate for the object. 

The latest-type object, HD 91782 B, resulted in a temperature much too high for its esti- 
mated spectral type of M9 according to all temperature scales and an unexpectedly low value 
for surface gravity, 3.02 dex, at the lower boundary of the surface gravity range of the models. 
This low surface gravity would imply an <1 Myr age, which is highly unlikely. The low grav- 
ity could a lso imply a background giant, b ut the star has been shown to have common proper 
motion by lMetchev fc Hillenbrand! (120091 ). It should be noted that for the latest type objects 
in particular, matching the near-infrared spectrum to a template spectrum is not a deter- 
mination of spectral type in the strictest sense and that different wavelength regimes and 
even spectral resolutions could imply different physical parameters, as is the case for numer- 
ous late-M a nd L dwarf s In pa rticular, the M8.5 7 dwarf 2MASS J06085283-2753583, was 
described by lRice et al.l (l2010ai ) as 2200 K and log(g)=3.20 dex [cgs] with a moderate resolu- 
tion (R ~2,000) J-band spectrum, but high-resolution spectra (R ~20,000) produced 2529 K 
and log(g)=3.98 dex [cgs], much more consistent with the age implied by confirmed member- 
ship in the Pictoris Moving Group and similarities to the well characterized TW Hydrae 
member 2MASS J12073347-3932540 jRice et aliboiObh . 



5.3. Case study: HD 129333 B 

Figure [7] shows the best-fit model spectrum plotted with the observed P1640 spectrum 
of HD 129333 B. The best fit parameters for HD 129333 B are T e/// =3506±211 K and 
log(g)=4.90±0.46dex [cgs]. The temperature is consistent with a M2.5±1V object according 
to spect ral type-effective tem perature conversions for either dwarf and intermediate gravity 



objects (ILuhman et al.lll999l ). along with the surface gravity matches the evolutionary tracks 
for a ~100 Myr, ~0.4 M Q object. These best fit parameters are plotted at the solid red line 
in Figure [7]). Also plotted are the 1-cr surface gravity range (red dashed lines) at the best bit 
effective temperature and the l-a effective temperature range at the best fit surface gravity 
(blue dotted lines). All spectra are normalized to the observations over the entire wavelength 
range. At these relatively high temperatures, the shape of the low resolution near-infrared 
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spectrum changes little even with ±200 K changes in temperature, but surprisingly the J- 
band slope seems somewhat sensitive to surface gravity, although all of the plotted models 
essentially fall within the uncertainties of the observed spectrum. 

The parameter distribution provided by the MCMC routine are more informative about 
the potential of this analysis. Figure [8] shows density contours of the accepted parameters 
from 9xl0 5 links in the MCMC. The highest contours in effective temperature are narrow 
and symmetric, providing a reliable (but large) 1-er uncertainty of 211 K. Although this is 
perhaps unpleasantly large, it is much smaller than the initial range of possible effective 
temperatures, 1400 to 4500 K, and did not require assumptions about the properties of the 
companion or the primary. The histogram for surface gravity, on the other hand, allows a 
large range of values, from ~4.1 to 6.0 dex, the upper boundary of the model grid. The 
success of low-resolution near-infrared spectra at producing reliable temperature results, 
but not fairing so well in gravity is also shown by using simulated P1640 spectra created 
using templates from the IRTF spectral library (Rice et al., in prep.). Thus suggests that 
constraining the age of the system via detailed studies of the primary star will be necessary 
in order to reliably characterize directly-detected low-mass companions. 



6. Null Results 



In addition to the five stars discussed above, we also observed 14 other stars that were 
part of the FEPS sample. These stars were observed in order to detect any additional 
companions than those previously known. We detected the known companion to HD 108944 
at the edge of our field of view with only a part of the PSF was visible; we could not extract 
a spectra or measure accurate astrometry. For the rest of the stars, we did not detect any 
candidates in the our 3'.'84 field of view. The companions of the known binaries are outside 
of our field of view. Table H] lists the HD number of the stars, the WDS number (indicating 
a known binary) and the Besselian date of the observations. In several cases, the stars were 



observed on multiple dates and we 



suppression software (ICrepp et al. 
same as that shown in 



rave listed all the dates. After data reduction with speckle 



Crepp et all (120 1 If ). 



20111). a ll the stars had contrast curves approximately the 



7. Summary 



We studied five binaries with late type companions. For each companion we present 
the astrometry from several epochs as well as an estimate of the spectral type derived from 
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Fig. 7.— Observed P1640 spectrum of HD 129333 B (black error bars) and synthetic 
PHOENIX spectrum (solid red) with best-fit parameters as determined by the spectral fit- 
ting procedure. The dashed lines (red in online version) shows two model spectra with 
T e fj=3506 K and surface gravities of 4.44 and 5.36 dex [cgs]. The dotted lines (blue in 
online version) shows two model spectra with surface gravities of 4.90 dex [cgs] and T e ff of 
3295 and 3717 K. The broken lines essentially demonstrate the 1-cr uncertainties listed in 
Table El 



matching measured spectra to template spectra. Most of the binaries have periods of the 
order of centuries, but HD 129333 already has a computed orbit and additional astrometry 
is needed to refine this orbit. This is particularly important, as the orbit is a combined 
RV and visual orbit, allowing computation of the individual masses of the stars. Additional 
observations will improve the mass estimates. HD 135363 also appears to have a period 
on the order of several decades and additional observations will allow an eventual orbit 
computation. These results increase our understanding of these particular systems, and also 
help to validate the P1640 data collection and reduction processes. We have also gained 
valuable experience in estimating the physical parameters of low mass companions using 
synthetic spectra calculated using atmosphere models. While potentially powerful, this tool 
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Fig. 8.— Density contours of Markov Chain Monte Carlo results for HD 129333 B in 
effective temperature (x-axis) and surface gravity (y-axis). 

is subject to systematic uncertainties both from the observed spectra (i.e., poor telluric 
correction or unreliable points at the ends of the spectrum) as well as from the models and 
intrinsic sensitivity of the spectra to the physical parameters (i.e,. lack of precision in best-fit 
surface gravity at low spectral resolution). 

In the future, we will collect multiple calibration observations in order to further reduce 
the spectral errors. In addition, during there have been advances in data reduction pipeline 
which promise to further reduce the noise of the extracted data. While it is being developed 
for the Phase II data, we will apply the new pipeline to the Phase I data. For Phase II data, 
improvements in the instrument have improved the optical throughput of the instrument, 
additionally decreasing the effect of noise. 
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Table 1. Measured Astrometry 



WDS 


HD 


Epoch 


on 




P(") 




Instrument 


09035+3750 


77407 


2002.0850 


353.364 ± 0.043 


1.6589 


± 


0. 


0036 


PHARO 






2003.0351 


353.796 ± 0.114 


1.6685 


± 


0. 


0032 


PHARO 






2003.3547 


353.890 ± 0.030 


1.6723 


± 


0. 


0029 


PHARO 






2009.2083 


355.843 ± 0.159 


1.6875 


± 


0. 


.0071 


PI 640 


10368+4743 


91782 


2002.1669 


33.572 


± 


0.455 


1.0022 


± 


0. 


.0084 


PHARO 






2003.0294 


33.551 


± 


0.265 


1.0189 


± 


0. 


.0028 


PHARO 






2003.3544 


33.833 


± 


0.097 


1.0250 


± 


0. 


0021 


PHARO 






2004.0969 


32.451 


± 


0.107 


1.0246 


± 


0. 


0019 


PHARO 






2004.4280 


32.532 


± 


0.104 


1.0322 


± 


0. 


0007 


NIRC2 






2004.4305 


32.639 


± 


0.106 


1.0243 


± 


0. 


0010 


NIRC2 






2009.2056 


30.5 


± 


2.1 


1.01 


± 


0. 


.08 


PI 640 


12547+2206 


112196 


2002.0876 


55.518 


± 


0.094 


1.5008 


± 


0. 


0013 


PHARO 






2003.0350 


54.839 


± 


0.095 


1.5114 


± 


0. 


.0012 


PHARO 






2003.3572 


54.710 


± 


0.086 


1.5157 


± 


0. 


.0016 


PHARO 






zUU4.Uy ( A 


54.351 


± 


0.085 


1.5241 


± 


0. 


UUzo 


TDTT A T? C\ 

r tlAKU 






2004.4851 


54.257 


± 


0.097 


1.5306 


± 


0. 


0009 


PHARO 






2009.2086 


50.62 


± 


0.26 


1.620 


± 


0. 


.008 


PI 640 


14390+6417 


129333 


2003.0296 


172.766 ± 0.114 


0.7172 


± 


0. 


0093 


PHARO 






2003.0323 


172.783 ± 0.102 


0.7286 


± 


0. 


0019 


PHARO 






2003.3629 


173.010 ± 0.318 


0.7361 


± 


0. 


0025 


PHARO 






2004.0974 


172.522 ± 0.178 


0.7415 


± 


0. 


0032 


PHARO 






2009.2007 


170.00 ± 0.55 


0.764 


± 


0. 


009 


PI 640 


15079+7612 


135363 


2002.0877 


121.352 ± 0.463 


0.2506 


± 


0. 


.0027 


PHARO 






2003.0378 


123.588 


- ± 0.354 


0.2564 


± 


0. 


.0024 


PHARO 






2004.4882 


127.384 ± 0.124 


0.2846 


± 


0. 


0007 


PHARO 






2009.2062 


133.14 ± 0.83 


0.363 


± 


0. 


.013 


PI 640 
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Table 2. The O - C Residuals for HD 129333 



UT 


A9 (°) 


Ap (") 


2003.0296 


-0.11 


0.027 


2003.0323 


-0.11 


0.027 


2003.3629 


0.18 


0.039 


2004.0974 


-0.21 


0.034 


2009.2007 


-2.09 


0.015 



Table 3. Spectral Results 



Name 


Sp. Type 


T e // (K) 


Log(g) [cgs] 


HD 77407 B 


M3V-M6V 


3322±141 


4.08±0.76 


HD 91782 B 


M9V 


3246±63 


3.02±0.07 [min] 


HD 112196 B 


M2V-M3V 


3914±320 


5.30±0.44 


HD 129333 B 


M2V-M3.5V 


3506±211 


4.90±0.46 


HD 135363 Ab 


M2V-M4V 


3667±281 


5.90±0.21 [max] 
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Table 4. Stars with No Detected Companions 



HD 


WDS 


Date (UT) 


70516 


08243+4457 


2009 1973 


90905 


10297+0129 


2009 2056 


92855 


10440+4612 


2009 2058 


95188 




2009.2058 


104860 


12046+6620 


2009.2087 


107146 




2009.2005 


108944 


12310+3125 


2009.2060 


134319 


15058+6403 


2009.2008 






2009.2035 






2009.2062 


145229 




2009.4929 


151798 


16501-1223 


2009.4903 


152555 




2009.2051 






2009.4875 


172649 


18397+3800 


2009.4904 






2010.5145 


203030 


21190+2614 


2008.8083 


204277 




2009.4881 






2009.4907 



